REPORT  DOCUMENTATION  PAGE  AFRL-SR-AR-TR-04- 

ning  the 

data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimE'.'e  or  any  <  /  J  Csf  educing 

this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (O'  CJ  L  r  /  2202- 

4302,  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  foi  »  currently 

valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

20-02-2004  FINAL 

3.  DATES  COVERED  (From  -  To) 

03-2001  to  12-2003 

4.  TITLE  AND  SUBTITLE 

Type- II  Quantum  Algorithms  for  Solitons 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

F49620-0lfD422 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

George  Vahala 

21 

3040303  199 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Physics 

College  of  William  &  Mary  College  of  William  &  Mary 

Williamsburg,  VA  23187 

8.  PERFORMING  ORGANIZATION  Kcrv^, 
NUMBER 

FINAL 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR ,  Directorate  of 

Mathematics  and  Space  Sciences 
-  Computational  Math. 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

Dr.  Jon  Sjogren 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT  - - - - 

There  are  many  problems  of  scientific  interest  that  are  not  tractable  on  any  foreseeable 
classical  computer.  Quantum  computers  have  the  potential  of  exploiting  non-classical 
features  like  quantum  entanglement  and  phase  coherence  which  can  exponentially  speed-up  the 
computational  algorithm.  Quantum  algorithms  have  been  developed  to  study  the  evolution  of 
solitons  in  the  Korteveg-de  Vries  equation  and  the  Nonlinear  Schrodinger  equation.  These 
nonlinear  equations  have  known  exact  analytic  solutions  to  which  the  quantum  algorithm 

solutions  have  been  compared  -  with  excellent  agreement.  Vector  soliton  propagation  down 

biref ringment  media  have  also  been  considered  as  well  as  soliton  turbulence  and  the 

corresponding  power  spectra.  It  has  been  found  that  two  on-site  qubits  per  spatial  node  is 

sufficient  to  model  a  scalar  continuum  field.  The  particular  choice  and  sequence  of  unitary 
collision  and  streaming  opereators  dictate  the  final  continuum  form  of  the  partial 
differential  equation. 

15  SUBJECT  TERMS  - - - - - 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

George  Vahala 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

757/221-3528 

BEST  AVAILABLE  COPY 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.1 8 


AFOSR  FINAL  REPORT 


"TYPE-II  QUANTUM  ALGORITHMS  for  SOLITONS" 


February,  2004 


George  Vahala 
William  &  Mary 


During  the  period  of  this  grant,  we  have  been  working  on  the  development  of 
Quantum  Algorithms  for  nonlinear  physical  systems,  in  collaboration  with  Dr.  Jeff  Yepez 
(Hanscom  Field)  and  Dr.  Linda  Vahala  (Old  Dominion  University).  In  particular,  quantum 
algorithms  have  been  developed  for 

(a)  solitons 

“Quantum  Lattice  Gas  Representation  of  Some  Classical  Solitons” 

G.  Vahala,  J.  Yepez  and  L.  Vahala 
Phys.  Lett.  A310,  187-196  (2003) 

“Quantum  lattice  gas  representation  for  vector  solitons” 

G.  Vahala,  L.  Vahala,  and  J.  Yepez 
SPIE  Conf.  Proc.  5105,  273  -  281  (2003) 

“Inelastic  Vector  Soliton  Collisions:  A  Quantum  Lattice  Gas  Representation” 

G.  Vahala,  L.  Vahala  and  J.  Yepez 
Phil.  Trans..  Roy  Soc.  London  (to  be  published) 

“Quantum  Lattice  Representation  of  Dark  Solitons” 

G.  Vahala,  L.  Vahala,  and  J.  Yepez 
SPIE  Conf.  Proc,  submitted  (2004) 


(b)  ID  MHD-Burgers  equation 

“Lattice  Boltzmann  and  Quantum  Lattice  Gas  Representations  of  One-Dimensional 
Magnetohydrodynamic  Turbulence” 

L.  Vahala,  G.  Vahala  and  J.  Yepez 
Phys.  Lett.  A306,  227-234  (2003) 


We  have  predominantly  concentrated  on  soliton  research  since  exact  solutions  are 
known  for  KdV  and  both  the  scalar  and  vector  nonlinear  Schrodinger  equation  (NLS)  as 
these  will  provide  a  stringent  test  on  our  quantum  algorithms.  The  spatial  dimension  is 
discretized  into  a  set  of  spatial  nodes.  For  modeling  either  the  KdV  equation  or  the  scalar 
and  vector  NLS  each  scalar  field  component  we  require  2  qubits  at  each  lattice  site.  The  on¬ 
site  qubits  are  entangled  by  the  unitary  collision  operator  and  this  entanglement  is  spread 
throughout  the  system  by  unitary  streaming.  In  particular,  the  KdV  equation  is  modeled  by 
the  tensor  product  of  the  on-site  unitary  collision  matrix 
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with  the  collide-stream  algorithm 


I  If/(t  +  At))  =  SXC+  St2  C.  Sf  C  S2  C,S l  C  S2  C+.S ,  CSl C  I y/(t)) 
where  Sx  is  the  global  streaming  operator  on  qubit- 1  to  the  neighboring  lattice  site  while  the 
transpose  streaming  operator  on  qubit-2  is  S2.  C  is  the  tensor  product  of  the  on-site  collision 
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matrix  UKdV  and  (7  is  the  adjoint  operator.  In  the  continuum  limit,  the  resulting  partial 
differential  equation  is  the  KdV  equation 
dy/  dy/  1  dV 

dt  dx  2  dx 3 

after  the  phase  transformation  is  introduced  to  yield  the  nonlinear  “potential”  term 
y/  -»  exp[zTAt]^  ,  with  V=i~- 

The  specific  streaming  sequence  is  required  to  eliminate  the  standard  diffusive/dispersive 
d2  /dx2  term  and  thus  give  the  required  leading  order  linear  term  of  the  KdV-equation 
c?  /dx3. 

To  recover  the  scalar  NLS  we  now  entangle  the  on-site  two  qubits  by  the  unitary 
collision  operator 
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and  utilize  the  global  collide-stream  algorithm 

|y<*  +  A*))  =  St2C  S2C  .St2CS2C  .  StxC  SxC  .SfCSx  C  \y/(t)) 

The  accuracy  of  the  finite  quantum  difference  algorithm  becomes  second  order  by 
symmetrizing  the  collide-stream  algorithm  on  each  on-site  qubit.  The  nonlinear  potential 
term  is  recovered  by  the  phase  transformation 

y/  — »exp[*TA t\y/  ,  with  F=|^|2 
to  yield  the  cubic  NLS  in  the  continuum  limit 

4r+  l-r+MV  -  o 


These  algorithms  are  readily  extended  to  consider  the  coupling  of  the  polarizations  due  to  the 
birefringent  medium.  In  this  case  we  have  two  coupled  NLS  equations,  requiring  two 
qubits/node  for  each  polarization.  The  coupling  between  the  two  polarizations  is  achieved  by 
the  appropriate  coupling  phase  transformation  that  now  couples  all  four  on-site  qubits.  An 
interesting  exactly  soluble  example  is  the  inelastic  collision  of  vector  Manakov  solitons  in 
which  for  one  of  the  polarizations  a  soliton  is  destroyed.  It  reforms  following  a  second 
collision 


